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ABSTRACT 



We propose near-infrared colour selection criteria to extract Active Galactic Nuclei (AGNs) using the near-infrared Colour-Colour 
Diagram (CCD) and predict near-infrared colour evolution with respect to redshift. First, we cross-identify two AGN catalogues 
with the 2MASS Point Source Catalogue, and confirm both the loci of quasars/AGNs in the near-infrared CCD and redshift-colour 
relations. In the CCD, the loci of over 70 ~ 80% of AGNs can be distinguished from the stellar locus. To examine the colours of 
quasars, we simulate near-infrared colours using Hyperz code. Assuming a realistic quasar SED, we derive simulated near-infrared 
colours of quasars with redshift (up to z ~ 11). The simulated colours can reproduce not only the redshift-colour relations but also 
the loci of quasars/AGNs in the near-infrared CCD. We finally discuss the possibility of contamination by other types of objects. 
We compare the locus of AGNs with the other four types of objects (namely, microquasars, cataclysmic variables, low mass X-ray 
binaries, and massive young stellar objects) which have a radiation mechanism similar to that of AGNs. In the near-infrared CCD, 
each type of object is located at a position similar to the stellar locus. Accordingly, it is highly probable that the four types of objects 
can be distinguished on the basis of the locus in a near-infrared CCD. We additionally consider contamination by distant normal 
galaxies. The near-infrared colours of several types of galaxies are also simulated using the Hyperz code. Although galaxies with 
Z ~ 1 have near-infrared colours similar to those of AGNs, these galaxies are unlikely to be detected because they are very faint. In 
other words, few galaxies should contaminate the locus of AGNs in the near-infrared CCD. Consequently, we can extract reliable 
AGN candidates on the basis of the near-infrared CCD. 

Key words, galaxies: active - galaxies: quasars: general - catalogs 



1. Introduction 

Active Galactic Nuclei (AGNs) are tremendous energetic 
sources, where vast amounts of energy are generated by gravita- 
tional accretion around supermassive black hole. The radiation at 
nearly all wavelengths enables us to detect AGNs in multiwave- 
length observations. Hence, AGNs have been studied at vari- 
ous wavelengths. Past studies show that their Spectral Energy 
Distributions (SEDs) are roughly represented by a power-law 
(i.e., f v cc v~ a ), whilst normal galaxies produce an SED that 
peaks at ~ 1.6/vm as the composite blackbody spectra of the 
stellar population. Because the colours of an object provide us 
with rough but essential information about its spectrum, colours 
are important clues to identify AGNs from normal stars. 

Colour selection is an efficient technique to distinguish 
AGNs from normal stars and have played an important role to 
extract AGN candidates without spectral obser vation. A clas- 
sic method is known as the t/V-excess (UVX; ISandagdfl96l 
ISchmidt & Greenlll983tlBovle et alJll990l) . The UVX technique 
exploits the fact that quasars are relatively brighter than stars 
at shorter wavelength and therefore occupy a bluer locus in 
a CCD with respect to stars. In addition, many AGN candi- 
dates have been sele cted on the basis of colours in various 
wavelengt hs: optical dRichards et alj 120021) . optical and near- 
infrared dGlikman et alj 120071) . and mid-infrared dLacv et alj 
2004; IStern et al.1 120051) . These studies provide us with clues 
about the properties of AGNs. 

Target selection of high redshift quasars has also been per- 
formed using their col o urs, m ainly in optical wavelengths (e.g., 
iFan et alj l2000l |2001L 120031) . However, near-infrared proper- 
ties are required when we try to select targets such as higher 



redshift quasars, since the shift of the Lyman break to longer 
wavelengths makes observations difficult at optical wavelengths. 
Therefore, near-infrared selection should be useful technique to 
extract high-redshift quasars. 

In this paper, we present a study of the near-infrared colours 
of AGNs and demonstrate, by both observed and simulated 
colours, that the near-infrared colours can separate AGNs from 
normal stars. Additionally, we predict near-infrared colour evo- 
lution based on a Monte-Carlo simulation. In Sect. |2l we in- 
troduce the catalogues of AGNs which are used in order to 
investigate the observed colours. We confirm the near-infrared 
properties of spectroscopically confirmed AGNs on the basis of 
the near-infrared CCD and redshift-colour relations in Sect. [3] 
In Sect. |H we simu late the near-infrared colours using Hyperz 
code developed by Bolzonella et al. (2000) and demonstrate that 
AGNs reside in a distinct position in the near-infrared CCD. In 
Sect. [5] we consider the other probable objects which are ex- 
pected to have near-infrared colours similar to those of AGNs. 



2. Data 

We examine the near-infrared properties of quasars/AGNs us- 
ing 2MASS magnitudes. The samples of quasars/AGNs are ex- 
tracted from Sloan Digital Sky Survey Data Release 5 (SDSS- 
DR5) quasar catalog and Quasars and Active Galactic Nuclei 
(12th Ed. ) and these catalogues are briefly introduced below. 
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Fig. 1. (a) The d istribution of AGNs in the OA catalogue , (b) The distribution of quasars in t he SO catalogue. The ste llar locus 
( Bes sell & Bretill 19881) . the CTTS locus dMever et al .11 19971) . and the reddening vector taken from lRieke & Lebofskvl d 1985b are also 
shown. 



2.1. 2MASS 

The T wo Micron All Sky Survey (2MASS Q ISkrutskie etaD 
2006) is a project which observed 99.998% of the whole sky 
at the J (1.25 yum), H (1.65 /mi), Ks (2.16 /an) bands, at Mt. 
Hopkins, AZ (the Northern Hemisphere) and at CTIO, Chile 
(the Southern Hemisphere) between June 1997 and February 



1 2MASS web site (http://www.ipac.caltech.edu/2mass/) 



2001. The instruments are both highly automated 1.3-m tele- 
scopes equipped with three-channel cameras, each channel con- 
sisting of a 256 x 256 array of HgCdTe detectors. The 2MASS 
obtained 4 121 439 FITS images (pixel size ~ 2"0) with 7.8 s 
of integration time. The 10cr point-source detection levels are 
better than 15.8, 15.1, and 14.3 mag at J, H, and K s bands. The 
Point Source Catalogue (PSC) was produced using these images 
and catalogued 470 992 970 sources. In the 2MASS web site, 
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redshift 
range 


QA catalogue 






SQ catalogue 




Region I 


Region II 


total 


Region I 


Region II 


total 


< z < 1 


1 671 (27) 


4 480 (73) 


6 151 


222 (11) 


1 869 (89) 


2 091 


1 < z< 2 


238 (47) 


265 (53) 


503 


222 (47) 


249 (53) 


471 


2<z<3 


67 (25) 


197 (75) 


264 


38 (19) 


165 (81) 


203 


3 < z < 4 


7(16) 


36 (84) 


43 


9(18) 


41 (82) 


50 


4 < z < 5 


5(71) 


2(29) 


7 


1(50) 


1 (50) 


2 


total 


1 998 (29) 


4 970 (71) 


6 968 


500 (18) 


2 317 (82) 


2 817 



Table 1. The number of objects distributed in the Region I and II. Of 7 061 AGNs detected in the QA catalogue, 93 do not have a 
measured redshift. The values in parentheses represent percentages of quasars/AGNs residing in each region. 



the images and the PSC are open to the public and are easily 
available. 



2.2. SDSS-DR5 quasar catalog 

The Sloan Digital Sky Survey (SDSS) provides a photometri- 
cally and astrometrically calibrated digital imaging survey of n 
sr above Galactic latitu de 30° in five broad optical bands to a 
depth of g' ~ 23 mag dYork et al.l2 000). Many astronomical cat- 
alogues have been produced by this survey. 

The SDSS quasar catalog IV dSchneider etal.ll2007l here- 
after SQ) is the forth edition of the SDSS quasar catalog I 
dSchneider et all 120021). which is made from the SDSS Fifth 
Data Release dAdelman-McCarthv et alj 120071) . The SQ cata- 
logue consists of 77 429 quasars, the vast majority of which were 
discovered by the SDSS. The area covered by the catalogue is 
«: 5740 deg 2 . The quasar redshifts range from 0.08 to 5.41, with 
a median value of 1.48. The positional accuracy of each object 
is better than 0"2. 

2.3. Quasars and Active Galactic Nuclei (12th Ed.) 

The catalogue of Quasars and Active Galactic Nuclei (12th Ed.) 
dVeron-Cettv & Veror3l2006l hereafter QA) is the 12th edition 
of the catalogue of quasars first published in 1971 by De Veny 
et al.. The QA catalogue contains 85 221 quasars, 1 122 BL Lac 
objects and 21 737 active galaxies (including 9 628 Seyfert 1). 
this catalogue includes position and redshift as well as optical 
brightness (U, B, V) and 6 cm and 20 cm flux densities when 
available. The positional accuracy is better than 1"0. 



3. Near-infrared properties of AGNs 

3.1. Extraction of Near-infrared counterpart 

The sources in two of the above-mentioned AGN catalogues 
(SQ and QA) were cross-identified with 2MASS PSC, and we 
extracted a near-infrared counterpart of each source. As men- 
tioned in the previous section, the positional accuracies in both 
catalogues are better than 1". Therefore, we identified an near- 
infrared counterpart when a 2MASS source is located within 1" 
of a SQ/QA position. 

As a result of the extraction, we have derived 9 658 (SQ cat- 
alogue) and 14 078 (QA catalogue) near-infrared counterparts. 
For investigating the near-infrared properties using 2MASS 
magnitudes, we used only 2 817 (SQ) and 7 061 (QA) ob- 
jects where 2MASS photometric quality flags are better than B 
(signal-to-noise ratio (S/N) > 7<x). 



3.2. Colour-colour diagram 

Near-infrared (H - K$)-(J - H) CCD is a powerful tool to inves- 
tigate the property of celestial objects. We investigated the near- 
infrared properties of quasars/AGNs using near-infrared CCD. 

Figure [TJ shows the distributions of quasars/AGNs in a (H - 
Ks)-(J - H) CCD. In previous studies, the intrinsic loci of 
stars and Classical T T auri Stars (CTTS ) were well defined by 
iBessell & Brettl d!988l) and iMever etaf] d 19971) . Their loci are 
also shown in the CCD. IBessell & Brettl d 1988b and the Caltech 
(CIT) systems are transformed int o the 2MASS ph otometric sys- 
tem by the method intr oduced by Carpenter (2001). The redden- 
ing vector, taken from iRi eke & Lebofs kv ( 1985), is also shown 
in the diagram. Because the stellar and CTTS loci can only shift 
along the reddening vector, most of these types of stars funda- 
mentally should not be located in the region described by the 
following equations. 

(J-H) < 1.70(//-tf s )-0.118 (1) 
(J -H)<0.6l(H-K s ) + 0.50 (2) 

Equation (JTJ represents the lower limit line where normal stars 
can reside and Equation (f2j) represents the lower limit line where 
CTTS can reside. Both lines are also shown in Fig. QJ Below, 
we call the region enclosed by Equations (QJ) and (f2| "Region II" 
and all other regions "Region I". In Fig.QJ we can see that most 
of the quasar/AGNs are located in clearly different areas than 
the stellar loci. The distributions of the quasar/AGNs are on the 
right side of the stellar loci in the CCD, i.e., they have a (J - H) 
colour similar to that of normal stars but have a (H - K$) colour 
redder than that of normal stars. Table [TJ counts the number of 
objects in each region. It shows that 70% of AGNs and 80% 
of quasars are distributed in Region II. Hence, the near-infrared 
selection of quasars can be more effective than that of other types 
of AGN. Especially, ~ 90% of low redshift quasars with < 
Z < 1 reside in Region II, so these quasars are rarely missed. 
However, objects with 1 <z<2or4<z<5 tend to have a 
bluer colour in (H - K$) than objects with other redshift ranges, 
which is similar to the colour of normal stars. Therefore, some 
of these quasars/AGNs might be missed. The difference of the 
loci between quasars/AGNs and normal stars is probably due to 
the difference of the radiation mechanism because the dominant 
radiation of quasars/AGNs is not blackbody radiation. 

This colour prop erty is considered to be caused by K-excess 
dWarren et al.l2 000). They proposed a KX method where quasars 
with a [V - J) colour similar to that of stars would be redder in 
(J - K) colour. In other words, this KX method can separate 
quasars and stars on the basis of their colours. This technique 
has been used for selec t ing quasar candida tes (e.g., ISmail et ail 
2008: Uurek et al.ll2008h iNakos et al.ll2009l) . The present work is 
a variant of the original KX technique, using the (J - H) versus 
(H - K§) diagram. 
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Fig. 2. Colours versus redshift for SDSS quasars. The redshifts 
are taken from the SQ catalogue. The red solid lines show the 
average colour evolutions with respect to redshift. 



Ay 


J-K s 


J-H 


H-K s 





0.78 (0) 


0.94 (0) 


0.86 (0) 


1 


0.94 (0) 


0.69 (0) 


0.60 (0) 


2 


0.26 (17) 


0.29 (9) 


0.14 (84) 


3 


0.79 (0) 


0.82 (0) 


0.59 (0) 


4 


1.0 (0) 


1.0 (0) 


0.88 (0) 



Table 2. Results of a KS test between average colour evolution 
and simulated colour evolution. The decimal values represent 
KS distance between two data and the values in parentheses rep- 
resent significance level (percentage) for each KS test. 



3.3. Colours versus redshift 

In Fig. |2] we plot the SDSS quasars, in three colours versus red- 
shift together with average colour evolutions with respect to red- 
shift. The redshifts are taken from the SQ catalogue. 

Each colour undergoes only a small change or dispersion 
with redshift. This is probably due to a variety of spectral shapes 
and/or a variety of extinctions. In the near-infrared CCD, this 
small colour change causes a small difference of AGN locus. 

These properties can be reproduced by the simulation as 
mentioned below. 

4. Simulating the near-infrared colours of quasars 

In this section, we demonstrate that the locus of quasars is well 
separated from that of normal stars on the basis of a simulation 
using a realistic SED of quasars. 

In order to simulate the near-infrared colours of quasars, 
we performed a Mon te-Carlo simulation with Hyperz code 
dBolzonella et a l. 2000). The Hyperz code can calculate photo- 
metric redshift based on an inputted spectral template library, 
which finds the best fit SED by minimizing the x 1 derived from 
the comparison among the observed SED and expected SEDs. 



The reddening effects are taken into account according to a se- 
lected reddening law. Although this code is usually used for 
estimating photometric redshifts, we use it to derive the near- 
infrared colours at various redshifts. 

First, we made a magnitude list containing randomly gener- 
ated J, H, Ks magnitudes, ranging from 8 to 16 mag (roughly 
coincident with a reliable range of 2MASS magnitude) and pro- 
duced 100 000 data sets. These data sets were subjected to SED 
fittings usin g the Hyperz c ode. A realistic SED o f quasars was 
taken from Ipolretta et ail d2007l) (i.e.. QSOl in jPolletta et all 
d2007l) are used). According to lPolletta et all d2007l) . the SED of 
the QSOl is derived by combining the SDSS quasar composite 
spectrum and rest-frame IR data of a sample of 3 SDSS/SWIRE 
quasars dHatzim inaoglo uet alj|2005l) . We used the reddening 
law from ICalzetti et al] (2000), which is prepared by default in 
Hyperz code. Inputting the data sets into the Hyperz code, we 
derived photometric redshifts with the probabilities associated 
with the value of x 1 - We only used objects having probabilities 
of > 99%. 

Figure[3]shows the simulated colour evolutions with respect 
to redshift. The curves in each diagram represent the simulated 
colours with Ay = 0, 1, 2, 3, 4 (from bottom to top), respectively. 
To find the best fits for the average colour curves, we performed 
Kolmogorov-Smirnov (KS) tests between average colour curves 
and each simulated colour curve. Table|2]shows the result of the 
KS tests. In all three colours, the colour evolution with Ay = 2 is 
the best fit among five Ay values. In addition, the redshift-colour 
relations of SQ quasars can be roughly reproduced by simulated 
curves with < Ay < 3. A variety of extinctions probably gen- 
erate the dispersion of the colours. It should be noted that both 
(/ - H) and (J - K$) colours steeply increase over z ~ 9. This 
is due to shifting the Lyman break over the J-band wavelength. 
This property can be useful for extracting high-redshift quasars. 

In Fig. |U the simulated colours with Ay = 2 are shown in 
the (H - K S )-(J - H) CCD, tracked by redshift evolution. An 
important point is that the simulated position is well separated 
from the stellar locus, that is, it is consistent with the loci of 
quasars/ AGNs shown in Fig. Q] A variety of extinctions causes 
the dispersion of the simulated position and this can probably 
reproduce the dispersion of the loci of quasars/AGNs in Fig. Q] 
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Fig. 3. Simulated colours versus redshift. The curves represent the simulated colour evolutions with Ay = 0, 1, 2, 3, 4 (from bottom 
to top), respectively. The SDSS quasars (left panel) and the average colour evolution (right panel) shown in Fig. [2] are also plotted 
in the diagram. 





Microquasars 


CV 


LMXB 


MYSO 


I 


16 (84) 


245 (75) 


11 (73) 


27 (93) 


II 


3(16) 


82 (25) 


4(27) 


2(7) 


total 


19 


327 


15 


29 



Table 3. The number and percentage of objects distributed in 
each region. The values in parentheses represent percentage. 



It is also consistent with the fact that the quasars with < z < 
1 have relatively redder colour in (H - K§) compared with the 
quasars with 1 < z < 2. 

Although it is difficult to distinguish high-redshift quasars 
among z < 8, we can extract high-redshift quasar candidates 
with z > 8 on the basis of a (H - K^)-(J - H) diagram because 
the ( J - H) colour steeply increases overz ~ 8. 

5. Discussion 

5. 1 . Other probable objects 

Although the locus of AGNs in the near-infrared CCD is differ- 
ent from that of normal stars, other types of objects might be 
distributed in the locus with similar properties to AGNs. If a po- 
sition in the CCD depends on the radiation mechanism, other 
objects with radiation mechanism similar to AGNs are also ex- 
pected to be located at the same position. Below, we further ex- 
amine the loci of four types of objects which have non-thermal 
radiation or which are considered to be bright in both near- 
infrared and X-ray wavelengths: Microquasars, Cataclysmic 
Variables (CVs), Low Mass X-ray Binaries (LMXBs), and 
Massive Young Stellar Objects (MYSOs). 

Sample objects are extracted from three catalogues, namely 
Microquasar Candidates (Microquasars; iCombi et al.l [2008) 
Cataclysmi c Binaries, LMXBs, and related objects (CVs and 
LMXBs; iRitter & Kolbl l2003bl) . and Catalogue of massive 



young stellar objects (MYSOs; IChan et al.l 1 1996b . First, we 
cross-identified each catalogue with 2MASS PSC, an d ex- 
tracted the near-infrared counterparts. ICombi et al.l ([2008) have 
cross-identified their catalogue with the 2MASS catalogue by 
adopting a cross-identification of 4". T he positional accurac y 
in Ritter & Kolb catalogue are ~ 1 " (IRitter & Kolbl l2003a>. 
The objects in the MYSO catalogue were selected from the 
Infrared Astronomical Satellite (IRAS) PSC whose typical po- 
sition uncertaintie s are about 2" to 6" dBeichman et al.l 1988; 
lHelouetalJ [T988). Therefore, we set positional criteria for the 
cross-identification to < 2" (CV and LMXB catalogues) and < 
4" (Microquasar and MYSO catalogues). We used objects with 
a 2MASS photometric quality better than B (i.e., S/N > Icr). 
Using 2MASS magnitude, they were plotted in a (H-Ks)-(J-H) 
diagram. 

Figure [5] shows the CCD of each object. In every case, a few 
objects are distributed around the locus of the AGNs, although 
most of the objects are distributed around the stellar locus or 
reddened region from the stellar locus. Table [3] lists the number 
and percentage of objects distributed in each region. Although 
the ratios of CVs and LMXBs residing in Region II are rela- 
tively larger than the other two types of objects, it is not more 
than ~ 25%. In addition, few objects have (H - K$) ~ 0.8 in 
Region II though most quasars/AGNs have this colour (see Fig. 
[T). Accordingly, contamination by these four types of objects 
should be a small fraction. 

This means that the dominant radiation of the four objects 
should be thermal radiation. The AGNs also radiate thermal ra- 
diation, but it is a very small fraction compared with the non- 
thermal component produced by accretion around supermassive 
black holes. Therefore, AGNs should be well separated by these 
four objects using the near-infrared colours. 
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Fig. 5. The distribution of four types of objects: Microquasars (upper left), Cataclysmic variables (upper right), Low Mass X-ray 
Binaries (lower left), and Massive Young Stellar Objects (lower right). The stellar locus and the reddening vector are same as in Fig. 
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The boundary between Region I and II is also drawn in the dia- 
gram. 



5.2. Contamination by normal galaxies 

Distant galaxies that appear as point-like sources might contam- 
inate the AGN locus in the near-infrared CCD. We confirmed the 
locus of normal galaxies in the near-infrared CCD by perform- 
ing a Monte-Carlo simulation as in Sect. [4] The SEP templates 
we used are seven spiral galaxies in lPolletta et alj d2007l) : spirals 
range from early to late types (SO-Sd). 

Figure [6] shows the simulated intrinsic colours (i.e., Ay = 0) 
of the seven galaxies. Galaxies with < z < 0.8 have intrinsic 
colours similar to those of normal stars (i.e., they are in Region 
I). Galaxies with 1.4 < z < 3 are distributed around the red- 
dened region of either normal stars and/or CTTS. Therefore, 
they should not be mistaken for AGN candidates. On the other 
hand, simulated colours with z ~ 1 are located in Region II. A 
fraction of AGN in Region II are possibly mistaken for galaxies 
with z ~ 1 • However, galaxies at z ~ 1 should not have enough 
brightness to be detected with mid-scale telescopes. Even the 
bright est galaxy has no more than M ~ -23 mag at SDSS r- 
band dBlanton et al.ll200U iBaldrv et al.ll2004l) . If such a galaxy 
were located at z ~ 1, its apparent magnitude would be m > 20 
mag at J-band. In addition, the apparent magnitude would be 
even fainter because most galaxies have M > -23 mag and the 
apparent brightness suffers extinction. Accordingly, only large- 
scale telescopes can observe these galaxies. Hence, few galax- 
ies should contaminate the AGN locus in the near-infrared CCD 
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with respect to the data where limiting magnitude is below 20 
mag. 

6. Summary and Conclusion 

We confirmed the loci of catalogued quasars/AGNs in a (H - 
Ks)-{J - H) diagram, of which over 70 ~ 80% are clearly sep- 
arated from the stellar locus. In addition, we simulated the near- 
infrared colours of quasars on the basis of a Monte-Carlo sim- 
ulation with Hyperz code, and demonstrated that the simulated 
colours can reproduce both the redshift-colour relations and the 
locus of quasars in the near-infrared CCD. We also predicted the 
colour evolution with respect to redshift (up to z ~ 11). Finally, 
we discussed the possibility of contamination by other types 
of objects. The locus of AGNs is also different from those of 
the other four probable types of objects (namely, Microquasars, 
CVs, LMXBs, MYSOs) that are expected to be located at a sim- 
ilar locus. We also demonstrated by a Monte-Carlo simulation 
that normal galaxies are unlikely to contaminate the locus of 
AGN s in the nea r- infrar ed CCD. 

iHewett et alj d2006l) investigated near-infrared colours of 
quasars using an artificial SED, but we proposed near-infrared 
colour selection criteria for extracting AGNs and studied both 
observed and simulated colours with quantitative discussions. 
An important point is that our selection criteria require only 
near-i nfrared photometric data, al though some previous studies 
(e.g., iGlikman et al.ll2007l 12008) used colour selections based 
on colours between near-infrared and optical wavelengths. In 
other words, our selection criteria make extraction of candidates 
easier because only near-infrared colours are needed. This tech- 
nique should also be useful when we search for high-redshift 
quasars, since they become very faint in the optical wavelength 
due to the shift of Lyman break. 

This paper demonstrates that near-infrared colours can be 
useful to select AGN candidates. If an additional constraint 
is imposed, more reliable candidates can be extracted. When 
we use the near-infrared colour selection with an additional 
constraint for near-infrared catalogues containing sources dis- 
tributed in large area (e.g., 2MASS, DENIS, UKIDSS, and fu- 
ture surveys), a lot of AGN samples (possibly over ~10 000) 
are expected to be d e rived i n a regio n over ~10 000 deg 2 . 
Kouzuma & Yamaokal (l2009al) (see also iKouzuma & Yam aoka 
2009b]) cross-identified between the 2MASS and ROSAT cata- 
logues, and extracted AGN candidates in the entire sky using the 
near-infrared colour selection in this paper. These large number 
of samples may provide us with clues about such as the evolu- 
tion of AGNs and X-ray background. Additionally, in our simu- 
lation, quasars with z > 8 can be extracted on the basis of near- 
infrared colours. This property might be helpful to search for 
high-redshift quasars in the future. 
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